Fano resonances at light scattering by an obstacle 
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It is shown that elastic resonance scattering of light by a finite-size obstacle with weak dissipation 
is analogous to quantum scattering by a potential with quasi-discrete levels and exhibits Fano reso- 
nances. Localized plasmons (polaritons) , exited in the obstacle by the incident light, are equivalent 
to the quasi-discrete levels, while the radiative decay of these excitations plays exactly the same 
role as tunnelling from the quasi-discrete levels for the quantum problem. Mie scattering of light by 
a spherical particle and an exactly solvable discrete model with nonlocal coupling simulating wave 
scattering in systems with reduced spatial dimensionality are discussed as examples. 
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Introduction and motivation. Light scattering by an 
obstacle is one of the fundamentalproblems of electro- 
dynamics, see, e.g., monographs (HHf. Nowadays interest 
to this problem increases even more than ever owing to its 
numerous applications in subwavelength optics, informa- 
tion processing, nanotechnologics and related fields 
Recent studies of resonant scattering by small parti- 
cles with weak dissipation rates (see [H Q and ref- 
erences therein) have revealed new and unexpected fea- 
tures, namely the inverse hierarchy of optical resonances, 
a complicated near-field structure, unusual frequency and 
size dependencies etc., which allow to name such a scat- 
tering anomalous. However, despite the highlighted pe- 
culiarities, the physical nature of both normal resonant 
(Rayleigh) and anomalous scattering is the same and 
may be described briefly as follows. Incident light excites 
localized electromagnetic modes in a scattering particle 
(plasmons or polaritons) oscillating with the frequency 
of the incident wave to. The corresponding oscillations 
of polarization of the particle result in an emission of 
electromagnetic waves with the same frequency to which 
constitute the scattered light. Resonances in this picture 
correspond to the cases when to occurs close to (the real 
part of) the eigenfrequency of one of the localized modes. 

In the present Letter we reveal that this physical pic- 
ture is analogous to the well-known case of Fano reso- 
nances in quantum physics 0] • This analogy sheds a new 
light to the phenomenon. It allows to employ powerful 
methods developed in the theory of the Fano resonances 
(such as, e.g., the Fcshbach-Fano partitioning theory) to 
describe the resonant light scattering. It also explains 
easily certain features of the anomalous scattering and 
related problems, namely sharp changes in the scattering 
diagrams (from preferably forward to backward) upon 
small changes in u>, recently observed in numerical sim- 



ulations [8( and a typical asymmetric shape of the res- 
onance lines (see Fig. [T]) obtained in the present Letter 
by analysis of the exact Mie solution of the light scat- 
tering problem by a spherical particle [H, Q. We also 
introduce and study an exactly solvable model account- 
ing for all main features of the resonance light scattering 
by a particle and exhibiting distinctive Fano resonances. 
The model simulates wave scattering in an array of cou- 
pled waveguides 0, photonic crystals and related 
systems. Remarkable similarity in manifestation of the 
resonances in all these cases helps us revealing fundamen- 
tal links between these phenomena. 

It should be stressed that though the Fano resonances 
in light scattering have been discussed before, attention 
has been paid generally to microscopic aspects of inelastic 
scattering in a simple plane geometry (see, e.g., Ref. [ill]). 
In contrast to that, we focus on elastic light scattering 
by a finite-size particle when diffraction plays a funda- 
mental role in the scattering process. The scattering 
is accompanied by a two-step transformation: incident 
plane wave — > localized electromagnetic modes — > scat- 
tered light. On account of the second step of this trans- 
formation (radiative damping) the localized modes have 
a finite life-time, being actually quasilocalized even at 
zero dissipation rate. Note however, that as long as the 
scattering of a continuous wave is concerned, the losses 
of energy by the quasilocalized modes are exactly com- 
pensated by the gain from the incident wave. Then, the 
amplitude of the modes becomes time- independent. 

Another route for the incident light is just to bypass 
the scatterer, as well as to excite (nonresonantly) local- 
ized modes whose eigenfrequencies lie beyond the close 
vicinity of lo. Interference of the incident and re-emitted 
light generates a complicated near-field structure and 
may give rise either to strong enhancement (constructive 
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interference) or strong suppression (destructive interfer- 
ence) of the electromagnetic field. The analogy to Fano 
resonances of a quantum particle scattered by a potential 
with quasidiscrete levels becomes straightforward. The 
two possible routes for the incident light correspond to 
the resonant and direct (background) scattering, while 
the radiative decay of the (quasi)localized modes is iden- 
tical to tunnelling from the quasidiscrete level. 

Light scattering by a small spherical particle. Though 
the analogy mentioned is pretty general and valid for 
any particle size and shape, for the sake of simplicity 
in what follows we consider light scattering by a small 
spherical particle described by the exact Mie solution [l|, 
|2fl. According to the solution for a plane polarized wave 
propagating along z-axis with vector E parallel to x-axis 
the intensities of waves scattered in a given direction are 
described by the expressions 



5« = J.W cos 2 V - 
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where the subscripts indicate the corresponding polariza- 
tion (relative to the incident wave), and up to a certain 
common multiplier 
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Here 9 and tp are the polar and azimuthal angles of 
the spherical coordinate frame whose center coincides 
with that of the particle and z-axis with the propaga- 
tion direction of the incident wave; pj (cos#) stands 
for the associated Legcndrc polynomial, P, (cos#) = 

dP x (cos 9)/d(cos 9) , scattering amplitudes ai and hi may 
be written in the following form: 



i> } {q,e)+iGf> (q,e) 



(a) 



(4) 



and similar for hi with replacement F^ a \G^ a ' — > 
pw^Q(b). w j iere e (^j jg the re lative (with respect to 
the environment) dielectric permittivity of the particle, 
q = kR, k stands for the wavenumber of the incident 
wave in vacuum and R for the particle radius. As for 
p(a,b) ^ Q(a.b) ^ they are expressed in terms of the Bessel 
[ J/ + i/2 (C)] an d Neumann [-/V) +1 / 2 (C)] functions. The 
corresponding expressions are cumbersome and may be 
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FIG. 1: (Color online) An example of typical Fano resonance 
profiles at elastic light scattering. Forward (blue) and back- 
ward (red) scattering by a small spherical particle calculated 
according to the exact Mie solution in the vicinity of the 
quadrupole resonance; the profiles for J|| S ' and 1^ are iden- 
tical, cf. Eqs. (|6]), 0. The media and constants model the 
sketch of experiment on the anomalous scattering proposed 
in Ref. [12], namely a colloidal potassium particle with ra- 
dius R — 6.2 x 1CP 6 cm immersed in a crystal of KC1. The 
refractive index of the crystal n K ci = 1-5; the dielectric per- 
mittivity of the particle is approximated by the Drude for- 
mula: e — 1 — LOp/(u + iv)(jJ, where uj p — 5.77 x 10 15 s _1 and 
v = vf/R is determined by collisions of free electrons with 
the particle surface; vf ~ 10 s cm/s. 



found in Ref. 



The partial scattering cross-section 
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We remind briefly results of the analysis of the Mie 
solution (for more details see, e.g., [i| and references 
therein) . Optical resonances are defined by the condition 

Gf 1 {q, q) = 0. For a small particle (q <C 1) the condi- 
tion G^p (q, ei) — 0, regarded as an equation for ep , does 
not have any real solutions, while solutions of the equa- 
tion G*[ a) (q, eg) = have the form e e = -{l+e)/£+0(q 2 ). 
Through the dependence e(uj), these solutions define res- 
onant frequencies u>i for non-dissipating materials. For 
weakly dissipating materials the roots of the equation 
(q,Ui) — have small imaginary parts. In this case 
the resonant frequencies arc equal to the real parts of the 
corresponding roots. Next, at q -C 1 amplitudes hi are 
small relative to ag and may be dropped. Regarding F^ , 
in this limit = O [q 2i+l ) and it does not vanish any- 
where but at the trivial point e = 1; as for G^ a \ away 
from the close vicinity of the resonance frequencies it is 
of order of one. All together that result in sharp optical 
resonances for cr S ca 1 with a typical symmetric Lorenzian 
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(Brcit-Wigner) profile [13j. 

(si 

However, such a profile is not always the case for 1^ 
Let us first discuss the non-dissipative limit (Im e = 
0). Considering excitations of non- resonant modes as 
background scattering and bearing in mind the nar- 
rowness of the resonance lines in this case Q we ob- 
tain that the corresponding amplitudes ae ~ aio, where 
aio = O (q 2l+1 ) is a purely imaginary quantity equal to 

~iF^ a \q, e)/G { £ a \q, e). Then, the dipolc mode {£ = 1) al- 
ways plays the dominant role in the background scatter- 
ing. Accordingly, in the vicinity of the dipole resonance 
the intensity of the dipole scattering always remains large 
relative to the contribution of the other modes. As a re- 
sult interference of the resonant scattering with the back- 
ground gives rise to only small corrections to the resonant 
mode and the profiles remain Lorcnzian. 

For higher-order resonances with I > 1 , which are also 
well pronounced in the non-dissipative limit even for a 
small particle the situation is qualitatively differ- 
ent. In this cases the resonance scattering arises from 
the corresponding background, whose amplitude ago is 
much smaller than that for the background dipole mode 
(aio). As e approaches q the resonant amplitude rapidly 
increases and becomes much greater than aio, reaching 
a maximum at e = eg and decreasing back to a^o upon 
further change of e. Thus, in the vicinity of eg there 
are two points (one to the left, the other to the right 
from eg) where \ag\ = \cl\q\. Therefore, in addition to the 
constructive interference of the resonant mode with the 
background dipole scattering at e = q, we also find anni- 
hilation of the two modes at one of the equal-amplitudc- 
points (passing the resonance adds tt to the phase of the 
resonant mode, so only one of the two points satisfies the 
necessary phase condition for the destructive interference 
with the background mode). In this case an asymmetric 
Fano resonance profile is obtained. 

To illustrate this reasoning let us consider the vicinity 
of the quadrupole resonance (£ = 2). Taking G^ 
at u> = u>£, employing expansions of the functions in pow- 
ers of q [3] and additionally expanding G% in powers of 

Se = e — 62 [we remind that G 2 a \q, £2) = 0] one obtains 
the following approximate expressions for /M : 
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There are two characteristic scales in Eqs. ©, ([3]), de- 
scribing (i) a sharp resonant lineshape centered at 8e = 
with width Y = q 5 /6 and (ii) full suppression of the 
scattering which occurs at much greater scale: Sen = 
-q 2 cos 29/24 cos 9 and 5e± = -q 2 cos 9/24, but still in- 
side the resonance region (l^e^jj -C 1). Note also that 
there is change of sign of <Jeii t ± with variation of 9 shifting 




FIG. 2: (Color online) Schematic view of a discrete ID chain 
of sites with a defect described by Eq. (J8j . Solid and dashed 
lines correspond to the resonant and background scattering, 
respectively, see the text for more details. 



of the annihilation point from one side of the resonance 
peak to another. 

Finite dissipation and/or increase in g will broaden 
the resonances. It results in an overlap of resonance lines 
of different orders, which eventually may produce very 
complex profiles. Nevertheless the phenomenon remains 
qualitatively the same as long as the dissipation is weak. 
Profiles of forward and backward scattering, calculated 
for a potassium colloidal nanoparticle immersed in a KC1 
crystal, are presented in Fig.[T]as an example. The calcu- 
lations are performed for a realistic dependence e(w) fit- 
ting actual experimental data. Note, that in accordance 
with Eqs. ([6|), ([7|) the points of destructive interference 
for the forward (9 = 0) and backward (9 = tt) scatter- 
ing lie^on different sides of the corresponding resonant 
peaks 
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Now let us inspect one-dimensional resonant wave scat- 
tering in systems with reduced spatial dimensionality 
Numerous examples of such systems could be found, e.g., 
in book [l(| • Fano resonances exhibited by these systems 
have exactly the same nature as those discussed above. 
To make it sure we consider a localized point defect in 
a simple ID discrete chain with interaction between the 
nearest and next-to-nearest neighbor sites (the nonlocal 
coupling is a fundamental feature of the model). The 
model is described by the following set of equations (see 
also Fig.©: 

^ + E c ™+A+;=°; <2=woWn; c^ = c», (8) 

3=-2 

where the integer n labels the spatial sites. All sites of 
the chain but a defect situated at n = are identical. 
Accordingly, C" ±1 = 1, C™ ±2 — 7 for |n| > 2. As for the 
defect, C±i — /U, C± 2 — M7 an d w o stands for the shift 
of the defect intrinsic frequency with respect to that for 
the other sites of the chain. 

The solvability condition for a linear travelling wave 
solution [<f> ~ exp(iuj q t — iqn)] in the defect-free region of 
Eq. ijHJ) results in the following dispersion relation: 



2(cosg + 7C0s2q), 



(9) 



where g is a (quasi)wave vector. To keep the one-to-one 
correspondence between u> q and q from the non-trivial 
range < q < tt we restrict I7I < 1/4. 

The defect-free region allows for another type of solu- 
tions whose permitted frequencies cover the entire band 



4 



of those for Eq. ©: -2(1 - 7) < uj q < 2(1 + 7). These 
solutions are obtained by assigning complex values to q in 
Eq. @q = 7r±zp; p > and correspond to exponentially 
decaying in space localized states. 

To study wave scattering by the defect in our model 
we consider an incident wave with amplitude equals one 
advancing from the left [l7| . Then, the corresponding 
boundary conditions read: 
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(10) 



where n and pi, ti and pi arc the transmission and 
reflection amplitudes of the propagating and evanescent 
waves from both sides of the defect, respectively. The 
transmission coefficient T is related to T\ by the expres- 
sion T = \ Tl \ 2 . 

The linear problem Eqs. (|5|)- (fTU|) is exactly solvable, 
but the solution is extremely cumbersome. It is remark- 
able, however, that the expression for n has a structure 
identical to that of Eq. ([4]), namely 



n = 



(11) 



where f ( Tl ) and are real quantities. Eq. (jTTJ) is 

simplified drastically at small 7, p and luq under the 
additional restriction of weak coupling between the defect 
and the chain p 2 <C \ll>qj\. In this case 
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where Suj = u> q — lu . The corresponding T{5u>) = 1 at 
Slu = 0, i.e., when ui q coincides with the intrinsic defect 
frequency and vanishes at 5u> = —p 2 /^, cf. the discussed 
above resonance scattering of light. 

To make sure that the obtained properties of the trans- 
mission coefficient are related to the resonant excitation 
of the defect, we calculate <f>Q. In the same approximation 
(small 7, p, ujq and p 2 <C |wo7|) we find: 



/' 



p 2 + 7(50-" — i5u)/2 



(13) 



Thus, \4>o\ 2 exhibits a typical sharp Lorenzian profile cen- 
tered about Suj w — 47/i 2 . 

More extensive analysis of Eqs. ([8l)- (fTU|) free from the 
imposed restrictions for 7, p and u>o does not add any 
qualitative difference to the results presented here, ex- 
cept for the fact, that the locations of resonant trans- 
mission and reflection may be separated strongly in the 
frequency space. Note, though the next-to-nearest sites 
interaction is crucial to the model, its generalization to 
more extended interactions (provided the coupling con- 
stants become smaller with increase in distance between 
the interacting sites) gives rise just to small quantitative 
corrections to the model characteristics. 



Conclusions. A deep connection between the light 
scattering by a finite obstacle and Fano resonances in 
quantum physics is revealed and illustrated by the anal- 
ysis of the exact Mie solutions and a one-dimenssion 
model. The concept developed in the present Letter may 
be generalized to more complicated problems, e.g., the 
scattering by a system of particles including coherent 
backscattering and weak localization phenomena, scat- 
tering in waveguides, etc. In each particular system, the 
specific resonant states arc different but all of them act 
in the same way being described by the similar physics. 
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